TRG-031-FR 


00 

Submitted  by: 

o 

TRG,  A  Division  of 

Control  Data  Corporation 

Route  110 

Melville,  New  York  11746 

AMPHIBIOUS  VEHICLE  STUDIES 
Under  Contract  No.  Nonr-4650(00) 


Distribution  of  this  document  is  unlimited. 


Submitted  to: 

Office  of  Naval  Research 
Department  of  Navy 
Code  438 

Washington,  D.C.  203G0 


„pc 

\>  V 


Reproduction  in  whole  or  in  part  is  permitted  for  any 
purpose  of  the  United  States  Government. 


June  1966 


TRG-031-FR 


Submitted  by: 

TRG,  A  Division  of 
Control  Data  Corporation 
Route  110 

Melville,  New  York  11746 


AMPHIBIOUS  VEHICLE  STUDIES 
Under  Contract  No.  Nonr-4650(00) 


Distribution  of  this  document  is  unlimited. 


Submitted  to: 

Office  of  Naval  Research 
Department  of  Navy 
Code  438 

Washington,  D.C.  20360 


Reproduction  in  whole  or  in  part  is  permitted  for  any 
purpose  of  the  United  States  Government. 


June  1966 


TABLE  OF  CONTENTS 


No. 


List  of  Figures . * . .  ii 

Acknowledgement  . . . . . . .  iii 

1.  Introduction . . . .  1 

2.  Analysis  of  LVTP5  Model  Test  Data .  -  2 

3.  Resistance  Calculations  . . . *....  8 

4.  Streamline  Calculations  for  LVTP5  Model . .  i . . . .  16 

4.1  Preliminary  Discussion .  16 

4.2  Normal  Rectangular  Source  Sheet  . .  20 

4.3  Simple  LVTP5  . . . . . .  20 

4.4  TRG  LVTP5  Model . . . . .  23 

5 .  Douglas  Resistance  Calculations  . .  31 

6.  The  Effect  of  Hull  Modifications  on  Resistance  ....  34 

7.  Conclusions  and  Future  Work  . . . .  33 

8 .  References  . . . . . . . .  39 


i 


SUBSIDIARY  OF  CONTROL 


DATA  CORPORATION 


LIST  OF  FIGURES 


Fig.  2.1  -  Total  Resistance  of  Three  LVTP5  Models  . . 

Fig.  2.2  -  Total  Resistance  Coefficients  of  Three  LVTP5 

Models  and  Various  Blocks  . . 

Fig.  2.3  -  wave  Resistance  Coefficients  of  Three  LVTP5 

Models  . . . . 


Fig.  3.1a  -  Resistance  Coefficients  of  Blocks  as  a  Function 

of  B/E,  for  Various  v/VgL,  With  L/H  Fixed —  r  9 

Fig.  3.1b  -  100  R^/jpv^L^  as  a  Function  of  B/H,  for  Various 

vA'gL,  with  L/H  Fixed  . .  10 

Fig.  3.2  -  Various  Measured  Total  Resistances  and  Calculated 

Wave  Resistances  for  LVTP5  Models  . . .  13 


Fig.  3.3  -  Various  Measured  and  Calculated  Wave  Resistance 

Coefficients  for  LVTP5  Models . .  14 


Fig.  4.1  -  Possible  Flows  Generated  by  a  Constant  Source 
Density  on  a  Rectangle  Normal  to  a  Uniform 
Flow . .  19 


Fig.  4.2  -  yz- Projections  of  Streamlines  Generated  by  a 

Constant  Source  Density  on  a  Rectangle  Normal 

to  a  Uniform  Flow . . .  21 


Fig.  4.3  -  Vertical  Sections  of  Body  Generated  by  a 

Constant  Source  Density  on  a  Rectangle  Normal 

to  a  Uniform  Flow .  22 


Fig.  4.4a  -  LVTP5  Body  Surface  used  for  Streamline  Calcula¬ 
tions  ,  Body  Plan  . .  24 

Fig.  4.4b  -  LVTP5  Body  Surface  used  for  Streamline  Calcula¬ 
tions  ,  Profile  . . . . . .  25 

Fig.  4.5  -  xy-Projections  of  Various  Calculated  Streamlines 

of  Flow  Past  LVTP5  Model  . . .  26 

Fig.  4.6  -  xz-Projections  of  Various  Calculated  Streamlines 

of  Flow  Past  LVTP5  Model  .  27 

Fig.  5.1  -  Various  Measured  Total  and  Calculated  Wave 

Resistances  for  LVTP5  Models  . -  32 

Fig,  7.1  -  Nominal  Propulsive  Efficiencies  of  Track- 

Propelled  X- 100  Model,  1250  lbs.  and  938  lbs. 

Model  Displacement,  Level  Trim .  37 


1ft/, 


SUBSIDIARY 


CONTROL  DATA  CORPORATION 


ACKNOWLEDGEMENT 


The  author  is  indebted  to  Mr.  A. I.  Raff,  Mr.  J.G.  Koelbel 
and  Mr.  G.  Sammis  of  TRG  for  their  assistance  in  preparing  this 
report,  and  to  Mr-  H.  Dugoff  of  Davidson  Laboratory  for  supplying 
resistance  data  on  Davidson  Laboratory's  model. 


1. 


1,  INTRODUCTION 

TRG's  program  calls  for  theoretical  and  experimental 
work  on  resistance  and  resistance  reduction,  and  on  improved  flow 
to  the  grousers.  We  began  by  examining  the  data  in  [2],  including 
the  components  of  resistance  ana  the  dependence  of  resistance  on 
test  conditions.  Also,  in  order  to  establish  a  connection  with 
test  data  other  than  that  on  LVTP5  models  we  correlated  that 
model  test  data  with  the  extensive  tests  of  Hay  [l]  on 
parallelopipeds .  This  work  is  described  in  Section  2.  We  then 
attempted  to  find  a  mathematical  model  (source  or  dipole  distri¬ 
bution)  of  the  LVTP5,  such  that  the  calculated  resistance  and 
body  shape  would  agree  sufficiently  well  with  LVTP5  data  to  serve 
as  a  basis  for  studying  shape  modifications  and  their  effect  on 
resistance  and  flow.  This  work  is  described  in  Sections  3  and  4. 
We  have  begun  to  investigate  bow  shape  modifications.  We  are 
aiming  at  the  design  of  compatible  devices  for  bow  modification, 
which  may  be  deployable,  inflatable  or  useable  for  fuel  or  other 
fluid  storage,  and  which  perform  their  hydrodynamic  tasks  without 
degrading  vehicle  functions.  This  work  is  described  in  Section  6. 
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2.  ANALYSIS  OF  LVTP5  MODEL  TEST  DATA 

Figure  2.1  shows  1/4- scale  model  resistance  in  pounds 
versus  model  speed,  full-scale  speed  and  Froude  number  for: 

a)  LVTP5  model  X-100,  1250  lbs.  displacement 
(80,000  lbs.  full  scale)  level  trim,  tracks 
run  at  zero  slip,  as  reported  in  [2].  The 
curve  is  slightly  faired  from  test  points. 

b)  LVTP5  model  built  and  tested  by  Davidson  Labora¬ 
tory,  full  scale  displacement  82,500  lbs.,  embody¬ 
ing  tread  envelope  but  no  treads  or  grousers. 

Model  built  to  1/12- scale,  data  expanded  to 
1/4-scale  by  Froude  scaling  neglecting  fric¬ 
tional  resistance.  Data  supplied  by  Davidson 
Laboratory. 

c)  LVTP5  model  X-100,  as  in  a)  above,  with  tracks 
stationary. 

We  use  Froude  scaling  neglecting  frictional  resistance  because 
frictional  resistance  is  estimated  to  be  less  than  10%  in  all 
relevant  cases  while  resistance  is  not  uniquely  defined  to  within 
10%  because  of  varying  test  conditions.  The  difference  between 
the  three  curves  in  Figure  2.1  is  believed  to  be  due  mainly  to  the 
different  treatment  of  the  tracks.  Scale  effect  and  minor  model 
shape  differences  are  believed  to  be  unimportant.  The  difference 
between  the  curves  is  substantial.  Tests  at  other  displacements, 
reported  in  [2],  also  show  a  substantial  variation  of  towing 
resistance  with  track  slip. 
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Figure  2.2  shows  drag  coefficient  =  R/^pv  A,  where 
A  is  static  submerged  frobtal  area,  for 

a)  LVTP5  model  X-100  as  described  in  a)  above, 

b)  LVTP5  model  built  by  Davidson  Laboratory,  as 
described  in  b)  above, 

c)  LVTP5  model  X-100  as  described  in  c)  above, 

d)  a  fully  submerged  flat  plate  whose  aspect  ratio 
is  that  of  the  LVTP5  (single  or  double  model,  being 
insensitive  to  aspect  ratio  in  this  range) , 

e)  a  surface- piercing  flat  plate  whose  static 
submerged  aspect  ratio  is  that  of  model  X-100  (B/H  =  2.40), 

f)  an  unrounded  block,  interpolated  and  Froude 
scaled  from  the  report  of  Hay  [lj,  B/H  =  2.40,  L/B  =  3. 

g)  a  block  having  rounded  forefoot  (forefoot 

radius/static  draft  =  .6  ),  interpolated  and  Froude- 

scaled  from  [1],  B/H  =  2.40,  L/B  =  5  > 

h)  as  f)  above  but  with  forefoot  radius/static 

draft  =2.4  ,  L/B  =  5 

With  reference  to  e) ,  f ) ,  g)  we  note  that  the  data 

depends  weakly  on  L/B  in  this  range  of  L/B.  Also  a) ,  b) ,  c)  use 
the  same  value  of  A. 

Examining  Figure  2.2  we  see 


a)  all  the  coefficients  are  speed-dependent, 

b)  rounding  a  block  in  various  ways  reduces  the 


resistance  by  a  factor  of  2  or  more. 


ZERO  SPEED 


We  tentatively  identify  the  quantity  CR  -  , 

the  speed- dependent  part  of  CRS  with  wave  resistance,  keeping  in 
mind  that  the  coefficients  of  friction  and  separation  drag  are 


independent  of  Reynolds  number  in  this  range  of  Reynolds  number 
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Figure  2.2 

Total  Resistance  Coefficients  of 
Three  LVTP5  Models  and  Various  Blocks 
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We  are  encouraged  to  make  this  identification  by  the  fact  that 
the  low-speed  asymptote  of  for  a  flat  plate  is  in  fact  equal 
to  Cp  in  the  absence  of  a  free  surface.  This  quantity  is  plotted 
in  Figure  2.3,  for  three  LVTP5  models.  We  see  that  the  result  * 
is  almost  identical  for  the  X-1G0  models  with  the  tracks  run  at 
zero  slip  and  held  stationary,  which  is  plausible.  The -Stevens 
results  are  much  lower.  We  are  unable  to  account  for  the 
difference  in  test  results  on  the  basis  of  differences  in  con¬ 
figuration  and  test  conditions.  Since  we  are  attempting  to 
predict  the  wave  resistance  theoretically  the  difference  is  vital. 
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3.  RESISTANCE  CALCULATIONS 

The  objective  of  these  calculations  is  to  help  determine 

the  validity  of  mathematical  models  (source  or  dipole  distributions) 

of  model  X-100,  as  regards  wave  resistance.  The  validity  of  the 

models  as  regards  the  shape  of  model  X-100  is  discussed  in 

Section  4.  Estimated  wave  resistance  coefficients  from'LVTP5 

model  tests  were  plotted  in  Figure  2.3 

We  are  looking  for  dipole  distributions  whose  calculated 

wave  resistance  approximates  the  measured  model  nrave  resistance. 

Since  the  theoretical  model  does  not  allow  for  tracks  and  since 

the  waves  made  by  this  blunt  vehicle  may  be  nonlinear  and  therefore 

not  calculated  correctly  we  do  not  seek  a  very  close  approximation. 

Before  proceeding  further  we  mention  the  result  of  an 

examination  of  the  effect  of  varying  beam  on  the  wave  resistance 

of  blocks,  keeping  draft,  length  and  speed  fixed.  In  thin-ship 

2 

theory  the  wave  resistance  varies  as  B  ,  while  for  ships  the 
measured  data  varies  like  Bn,  where  n  varies  but  is  generally 
between  1.5  and  2.  The  data  in  [1]  indicates  that  for  blocks 
having  L/H  =  5.72  the  resistance  is  roughly  proportional  to  B^'^ 
at  full  scale  speeds  above  4  MPH,  in  the  range  1.5  <  B/H  <  3.0  or 
.263  <  B/L  <  .526  (see  Figure  3.1).  This  is  a  measure  of  the 

extent  to  which  we  are  out  of  the  thin-ship  domain. 

Wave  resistance  calculations  were  made  for  the  LVTP5 
model  X-100,  having  the  following  characteristics: 

Displacement  (pgV)  1250  lbs.  (freshwater) 

Length  (L)  7.5* 

Beam  (B)  3.0' 


Draft  (H) 

6  =  H/L 

TIM  B/2L  =  .2 
IKU  /a  subsidiary 


1.25'  (to  hull  bottom) 
.167 


I 


2 


3 


4  MPH 


MODEL 


MPH  TULL  SCALE 


Figure  3.1A 

Resistance  Coefficients  of  Blocks  as  a  Function  of  B/H 
For  Various  v/VgL,  With  L/H  Fixed. 
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The  value  of  A/2L  ,  a  quantity  which  appears  in  the  calculations, 


was 

A/2L2  =  -ggfc--  =  . 142 
2HL  pg 


In  the  first  set  of  calculations  the  model  was  represented 

by  a  constant- strength  volume  distribution  of  horizontal 'dipoles. 

The  region  occupied  by  the  dipoles  was  a  box.  The  length  and 

draft  Hp  of  the  box  were  fixed  equal  to  the  model  length  and  mode.l 

draft  H,  while  the  beam  B^  of  the  box  was  varied  (the  subscript  D 

denotes  dipole) .  The  constant  strength  was  chosen  to  make  linearized 
volume  equal  to  model  X-100  volume. 


Case 

1: 

Box 

length 

=  model  length 

Box 

beam  = 

0  (Centerplane  distribution,  Michell  theory) 

Case 

2: 

Box 

length 

=  model  length 

Box 

beam  = 

1/2  model  beam,  (3  =  Bp/2L  =  .1 

Case 

3: 

Box 

length 

=  model  length 

Box 

beam  = 

model  beam,  (3  =  .2 

These  three  calculations  show  the  effect  of  spreading  the  dipoles 
laterally.  The  next  set  of  calculations  were  the  same  except  that 
the  stern  waves  were  eliminated  by  extending  the  model  to  infinity 
downstream.  A  possible  justification  for  this  is  that  the  X-100 
is  blunt  and  has  a  strong  wake,  so  that  stern  wave  generation  is 
small. 

Case  4:  Box  length  =  semi- infinite  downstream 
Box  beam  =  0 


Case  5: 


Box  length  =  semi- infinite  downstream 
Box  beam  --  1/2  model  beam,  (3  =  .1 
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Case  6:  Box  length  =  semi- infinite  downstream 
Box  beam  =  model  beam,  p  =  .2. 

The  resulting  wave  resistance  in  pounds  is  plotted  in  Figure  3.2, 
together  with  the  total  resistance  from  model  tests. 

The  calculated  results  for  Cases  1,  2,  3  are  much  too 
high,  and  also  display  oscillations  not  present  in  the  measured 
data.  The  calculated  results  do  decrease  as  the  distribution  is 
spread  laterally,  as  a  result  of  increased  wave  interference. 

The  calculated  results  for  Cases  4,  5,  6  show  that  by  eliminating 
stem  waves  we  have  reduced  the  mean  resistance  and  eliminated 
the  oscillations,  as  expected.  Also,  resistance  again  decreases 
as  the  distribution  is  spread  laterally.  Nevertheless,  the 
numerical  values  of  the  calculated  wave  resistance  are  still  too 
high,  even  for  3  =  .2.  Also,  comparing  against  total  measured 
resistance  the  trend  at  high  speed  appears  wrong. 

When  we  examine  calculated  wave  resistance  coefficients 

the  picture  darkens  somewhat  (Figure  3.3).  First  of  all 

there  are  the  inevitable  problems  at  very  low  speed.  For  3=0 

the  calculated  resistance  coefficients  approach  non-zero  constants 

as  the  Froude  number  goes  to  zero,  the  limiting  value  being  twice 

as  great  with  stern  waves  as  without.  If  3  >  0,  the  resistance 

coefficient  vanishes  like  as  f  ->  0.  Also,  the  calculated 
coefficients  are  too  high  and  have  the  wrong  trend  at  high  speed. 

The  six  dipole  distributions  tried  above  are  quite 

rough.  In  the  course  of  studying  the  flows  generated  by  them 
(see  Section  4)  it  was  found  that  it  is  not  a  good  idea  to  spread 
the  distribution  too  much  laterally  since  after  some  spreading 
the  flow  no  longer  generates  a  closed  body.  While  it  may  turn 
out  that  resistance  can  be  calculated  with  useful  accuracy  using 
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distributions  which  almost  or  approximately  generate  a  closed 
body  we  prefer  not  to  pursue  this  approach  till  we  know  more 
about  it,  especially  since  the  functional  form 
(f(x,y,z)  =  constant  in  a  box)  of  the  dipole  distribution  is 
itself  a  gross  approximation. 

In  Section  4  we  say  how  we  have  used  the  Douglas  Program 
[3]  to  calculate  the  flow  past  an  LVTP5  model.  Part  of  the 
Douglas  Program  output,  namely  the  quadrilateral  descriptions 
and  source  strengths,  can  be  used  as  inputs  to  our  wave  resistance 
program.  Wave  resistance  calculations  of  this  type  are  described 
in  Section  5. 


I L 
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4.  STREAMLINE  CALCULATIONS  FOR  LVTP5'  MODEL 
4.1  Preliminary  Discussion 

Among  the  simple  dipole  distributions  which  are 
candidates  for  representing  the  LVTP5  model  are  those  defined 
by  functions 

£j)(x3y3z)  =  fD(x)g(yJz) 

where 

^(x)  =  1  for  j x |  <  ^ 

-  0  elsewhere 


(4.1)  ' 

(A) 


or 


while 


f^x)  =1  for  x  £  ^ 
-  0  elsewhere 


(B) 


g(yA)  =  constant  if  y,z  in  the  region  Ag  in  the 
yz-plane 
=  0  elsewhere. 

(A)  represents  a  model  of  finite  length,  whose  stern  waves  are 
undiminished  by  the  wake.  (B)  represents  a  semi- infinite  body 
having  no  stem  waves.  The  product  form  (4.1)  implies  that  the 
body  will  tend  to  have  constant  vertical  sections  except  near 
its  ends.  The  form  of  this  vertical  section  depends  on  the 
region  Ag  and  the  function  g(y,z).  A  completely  equivalent 
description  is  via  a  source  distribution 

fs(x,y,z)  =  [6(x  -  |)-5(x  +-j)Jg(y,z)  (A) 

or 

fs(x,y,z)  -  6(x  -  \)  g(y,z) 
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where  g(y,z)  is  as  above  and  5(x)  is  the  Dirac  6-function. 

We  have  evaluated  resistance  using  both  (A)  and  (B) , 
and  discussed  the  results  in  Section  3.  We  now  consider  the  choice 
of  suitable  regions  Ag  and  functions  g(y,z)  so  as  to  approximate 
the  LVTP5  model.  Since  the  dipole  distributions  introduced  above 
are  too  simple  to  represent  the  details  of  the  LVTP5  model  we 
introduce  an  equivalent  box  by  requiring 


IT  «  L  =  7.5'  (E  denotes  equivalent) 

BE/HE  «  B/H  =2.40 

VE  =  V  =  LEBEHE  =  20  ft? 

BE  =  2.530 
HE  =  1.054. 

The  equivalent  box  has  the  same  length,  beam/draft  ratio  and 
displacement  as  the  LVTP5  model. 

A  box- like  dipole  distribution  is  appropriate  to  approxi¬ 
mate  a  box- like  body,  so  we  choose  the  region  Ag  to  be  a  rectangle. 
We  intend  to  compute  the  body  from  the  dipole  distribution,  using 
the  streamlines  of  the  zero-speed  flow  (i.e.  thi  flow  satisfying 
b0/^z  =  0  on  the  free  surface) .  In  effect  this  gives  us  the 
double  body  (body  plus  image  in  the  free  surface) .  Keeping  this 
in  mind  we  are  led  to  an  interesting  constraint  on  the  choice  of 
the  rectangle  Ag  on  which  a  constant-strength  source  distribution 
is  to  be  placed.  Considering  the  stemless  or  semi-infinite 
situtation,  the  time  rate  of  fluid  addition  is  47rOgAg  and  far 
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downstream  where  the  flow  speed  is  unity  the  area  A  required 

CO 

to  transmit  this  fluid  is 


A 

CO 


= 


There  are  two  flow  possibilities  depending  on  the  value  of 
as  =2^-  (A^/Ag)  ,  as  s*lown  in  Figure  4.1.  The  x-veloeity  u+  on 
the  upstream  face  of  is  given  by 


If  u+  =  2?r0g  <  1  then  at  every  point  in  space  the  total  x- velocity 
is  downstream,  and  there  is  no  stagnation  point  or  closed  body. 

If  u+  =  ^Fcrs  ^  ^  fchen  there  is  a  stagnation  point  and  a  closed 
body,  since  the  total  x-velocity  on  the  upstream  face  of  Ag  is 
upstream.  Hence 

2irOg  >  1 

or 


A  >  2AC 

oo  O 

is  a  necessary  and  sufficient  condition  for  the  generation  of  a 
closed  body  (ignoring  the  border  line  case  2iro^  -  1) .  This  means 
that  if  we  wish  to  approximate  the  flow  past  an  LVTP5  model 
(with  or  without  stern)  by  a  flow  of  the  type  constructed  above, 
and  with  the  approximate  flow  to  generate  a  closed  body,  then 
approximately  we  must  require  Ag  <  ^  A^,  where  for  A^  we  use  the 
LVTP5  cross-section  area  at  the  central  station.  Recalling  the 
resistance  calculations  of  the  previous  section,  3  =  .2  would 
mean  that  no  closed  body  is  generated,  while  3  =  .1  is  borderline. 
As  stated  elsewhere,  the  requirement  that  a  closed  body  be 


generated  may  not  be  mandatory. 
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u+=  2  xr  <rs  <  I,  NO  STAGNATION  POINT, 
NO  CLOSED  BODY. 


2  -  DIMENSIONAL  SKETCH  ILLUSTRATING 
3-  DIMENSIONAL  SITUATION 


Source  Density  on  a  Rectangle  Normal 
to  a  Uniform  Flow. 
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We  have  calculated  some  streamlines  of  flow  past  a 
normal  rectangular  source  sheet,  as  described  above,  which 
generated  a  closed  body.  They  are  described  below  in  Section  4.2. 
In  addition,  surface  source  distributions  were  obtained  for 
box- like  bodies  approximating,  in  two  steps,  the  LVTP5  vehicle 
configuration.  These  are  described  in  Sections  4.3  and  '4.4. 

4.2  Normal  Rectangular  Source  Sheet 

We  have  run  some  streamlines  which  confirm  that  if 
2ircrc  <  1  no  closed  body  is  generated.,  For  illustration  we  show 
some  streamlines  for  a  case  in  which  there  is  a  closed  body.  The 
bow  rectangle  Ag  is  at  x  =  0.5  and  is  defined  by  |yj  <  .113, 

|Z|  <  .0942,  Ag  =  0.01065.  We  have  Gg  =  0.177,  so  that 
A^  =  2.22  Ag.  The  stagnation  point  was  found  at  x  =  -f  0.5116278. 
Figure  4.2  shows  the  y-z  projections  of  streamlines  starting  at 
this  value  of  x  and  at  a  distance  of  0.001  from  the  x-axis  at 
various  indicated  angles  from  the  y-axis.  It  shows  that  all 
such  body  streamlines  are  tangent  to  the  x-y  plane  at  y  =  0.0. 
Figure  4.3  shows  sections  through  the  half-body  generated  at 
various  values  of  x. 

4.3  Simple  LVTP5 

A  simplified  body  close  in  shape  to  the  LVTP5  was  drawn 
up.  The  surface  was  divided  into  59  quadrilaterals  of  about 
uniform  area  (for  1/4  of  the  submerged  body)  and  the  constant 
source  strength  over  each  quadrilateral  was  found  (via  the 
Douglas  Program  [3])  which  satisfied  the  boundary  condition  at 


M/a  subsidiary  of  control 


l/  r\ 


ORPORATION 


STREAMLINES 
SHOWN  TO 
X  =  0.5 


ALL  STARTING  POINTS  AT  X  =  .5116278  AND  .001  FROM  X-AXIS. 
ANGLE  INDICATED  IS  THAT  OF  STARTING  POINT  FROM  Y  —  AXIS 


Figure  4.2 

yr-Projections  of  Streamlines  Generated 
by  a  Constant  Source  Density  on  a 
Rectangle  Normal  to  a  Uniform  Flow. 


Normal  to  a  Uniform  Flow. 
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one  point  on  each  quadrilateral.  Using  a  Douglas  Program  Post 
Processing  Program  developed  by  TRG,  and  TRG’s  Streamline  Program, 
various  streamlines  were  calculated  for  the  resulting  singularity 
distribution.  These  streamlines  were  poor  approximations  to  the 
true  flow  field  in  that  they  passed  through  the  body.  This  was 
caused  by  negative  sources  on  the  front  of  the  body.  These 
negative  sources  are  thought  to  be  the  erroneous  result  caused 
by  using  too  few  quadrilateral  elements  to  approximate  the  body 
surface.  These  results  indicated  the  nature  of  the  refinement 
necessary  to  obtain  reasonable  streamlines  from  approximate 
surface  singularity  distributions  determined  using  the  Douglas 
program  and  our  programs. 

4.4  TRG  LVTP5  Model 

A  surface  singularity  distribution  was  calculated  for 
the  LVTP5  hull  shape  (with  the  tracks  and  wheels  removed  and  the 
well  filled  in).  Figure  4.4  shows  the  assumed  body  surface 
which  was  divided  into  223  quadrilaterals  (over  1/4  of  the  sub¬ 
merged  body) .  The  quadrilaterals  were  smaller  in  the  bow  region 
than  towards  amidships.  The  Douglas  Program  was  applied  to  this 
input  body  surface  and  calculated  the  source  strength  for  each 
quadrilateral  element,  uniform  over  each  element,  which  satisfied 
the  boundary  condition  at  one  point  on  each  element.  This 
singularity  distribution  was  used  to  calculate  the  streamlines 
of  the  flow.  A  saddle  stagnation  point  (see  Reference  [4])  was 
found  on  the  x-axis  at  x  =  14.230".  This  type  of  stagnation  point 
is  believed  to  be  due  to  the  indented  V-bow.  Figures  4.5  and  4.6 
show  the  x-y  and  x-z  projection  respectively  of  various  calculated 
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Figure  4.4b 

LVTP5  Body  Surface  used  for 
Streamline  Calculations,  Profile 


Streamlines  of  Flow  Past  LVTP5  Model 


Streamlines  of  Flow  Past  LVTP5  Model 
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streamlines  in  relation  to  the  LVTP5  body  surface.  Here  it  is 
seen  that  streamlines  very  close  to  the  body,  such  as  the 
waterline  or  keel  line,  can  be  obtained  without  difficulty.  The 
skew  streamlines,  which  are  also  very  close  to  the  body  surface, 
indicate  the  predominant  lateral  flow  on  the  bow  ramp  of  the 
vehicle.  Since  the  streamlines  for  the  generated  singularity 
distribution  are  reasonable,  a  calculation  of  linearized  volume 
was  made  based  on  the  resulting  223  source  strengths.  The  result 
is  a  linearized  volume  of  447  cubic  inches  (per  1/4  of  submerged 
model).  The  actual  volume  is  395  cubic  inches.  This  ratio  of 
linearized  to  actual  volume  (447/395  =  1.13)  is  about  what  we 
expect,  based  on  the  known  relation  between  true  volume,  linearized 
volume  or  dipole  moment,  and  added  mass  in  the  forward  direction, 
and  is  a  partial  check  on  the  accuracy  of  the  calculation  of  the 
source  strengths.* 

A  sample  of  computer  input  data  and  a  sample  of  computer 
output,  for  streamline  calculations,  is  given  on  the  following  two 
pages . 


*  Previous  experience  with  shiplike  forms  generated  by  singularity 
distributions  also  indicates  that  a  ratio  of  1.13  is 
reasonable. 
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5.  DOUGLAS  RESISTANCE  CALCULATIONS 

The  application  of  the  Douglas  Program  to  the  idealized 
LVTP5  was  described  in  Section  4.  The  resulting  source  distri¬ 
bution  has  a  wave  resistance  which  can  be  calculated  using  the 
point  source  approximation,  in  which  each  quadrilateral  carrying 
constant  source  strength  is  replaced  by  a  point  source  (having 
the  same  total  strength)  at  its  centroid.  We  calculated  the 
wave  resistance  in  this  way,  without  stem  waves.  The  result 
is  plotted  in  Figure  5.1.  At  full-scale  speeds  above  about  3  MPH 
the  calculated  resistance  rises  rapidly  and  becomes  about  twice 
the  block  resistances  (f3  =  .0,.1,.2  in  Figure  5.1).  At  low 
speeds  this  calculated  resistance  falls  below  the  previous 
calculations  for  f3  =  .0,  .1.  The  high  values  can  be  explained 
on  the  basis  of  the  factor  of  2,  discussed  in  Section  4,  which 
relates  the  source  density  generating  a  stagnation  point  and 
the  source  density  generating  the  proper  volume.  We  expect  to 
find  the  higher  value  a  =  =  .159,  in  the  case  of  the  Douglas 

output,  near  the  bow  of  the  vehicle,  and  inspection  of  the 
printed  source  strengths  reveals  that  indeed  the  values  of  a  are 
in  the  range  . 13  -  .17.  According  to  this  reasoning  the  Douglas 
resistance  should  be  roughly  four  times  the  resistance  plotted 
for  D  =  .2  in  Figure  5.1,  and  it  is.  At  low  speeds  the  sources 
over  the  sides  and  bottom  of  the  vehicle,  which  are  present  in 
the  Douglas  calculation  but  not  in  the  block  calculations,  cannot 
be  neglected  and  the  reasoning  becomes  incomplete. 
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Figure  5.1 

Various  Measured  Total  and 
Calculated  Wave  Resistances  for  LVTP5  Models 
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On  tne  basis  of  a  comparison  made  on  a  Series  60  hull, 
for  which  wave  resistance  based  on  the  Douglas  Program  was 
calculated  exactly  and  using  the  point-source  approximation,  we 
believe  that  the  point- source  approximation  is  adequate  except 
at  the  lowest  speeds  of  interest.  In  view  of  the  unpromising 
nature  of  the  calculated  results  based  on  the  point-source 
approximation,  and  the  anticipated  cost  of  the  exact  calculations, 
we  decided  not  to  make  the  exact  calculations. 

It  appears  that  wave  resistance  calculations  based  on 
the  Douglas  Program  do  not  provide  useful  results  for  the 
idealized  LVTP5.  One  possible  explanation  is  that  wave  velocities 
are  not  used  in  the  Douglas  program  which  determines  the  source 
strengths,  so  that  the  source  strengths  are  not  right.  Another 
possible  explanation  is  that  a  theory  in  which  the  free  surface 
is  linearized,  even  with  the  proper  source  strength,  is  not  a 
useful  approximation  for  bodies  as  blunt  as  the  LVTP5. 

The  block  resistance  calculations  give  the  best 
resistance  predictions,  even  though  they  do  not  represent  the 
shape  of  the  hull  as  accurately  as  the  Douglas  source  distribution 
does . 
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6.  THE  EFFECT  OF  HULL  MODIFICATIONS  ON  RESISTANCE 

Using  the  block  model  of  the  LVTP5  resistance  we  have 
considered  bow  modifications  intended  to  reduce  resistance.  Not 
a  great  deal  was  done  in  this  direction  but  some  interesting 

results  were  obtained.  In  Figure  5.1  the  curve  — o — ro - o — - 

is  calculated  wave  resistance  without  stem  waves  for  the  block 
representation  of  the  LVTP5  (with  3  =  .2),  modified  by  a 
bulb.  The  bulb  is  generated  by  a  horizontal  line  of  horizontal 
dipoles,  the  dipole  strength  being  given  by  a  quadratic  function 
of  x  whose  coefficients  were  determined  by  digital  optimization 
at  a  Froude  number  f  =  .35.  We  see  that  the  calculated  wave 
resistance  is  indeed  reduced  at  f  =  .35,  and  also  at  other 
Froude  numbers.  The  dipole  distribution  representing  the  bulb 
is  rather  wild,  and  we  did  net  have  time  to  verify  that  a  closed 
body  is  in  fact  defined.  If  a  closed  body  is  not  defined  the 
next  step  would  be  to  determine  what  resistance  penalty  is 
incurred  by  restricting  the  optimization  to  realistic  dipole 
distributions.  The  large  calculated  resistance  reduction 
suggests  that  the  subject  is  worth  pursuing. 
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7.  CONCLUSIONS  AND  FUTURE  WORK 

The  practical  aims  of  the  project  are  to  produce  new 
configurations  which  will  have  increased  speed  in  the  water. 

We  have  begun  by  seeking  an  improved  theoretical  understanding  of 
the  measured  data  obtained  to  date.  Some  points  to  be  established 
are: 

1)  what  fraction  of  the  vehicle  drag  or  power 
consumption  is  attributable  to  wave  resistance, 

2)  what  fraction  of  the  vehicle  drag  or  power 
consumption  is  attributable  to  poor  flow  to  the 
foremost  grousers. 

Regarding  1)  we  have  found  (in  Section  2)  that  the  estimated  wave 
resistance  lies  between  15%  and  50%  of  the  total  resistance  at 
full-scale  speeds  in  the  neighborhood  of  7  MPH.  This  uncertainty 
is  large.  The  total  resistance  is  itself  uncertain  to  the  extent 
of  some  30%.  Regarding  2)  we  know  from  Figure  277  of  [2]  that 
at  model  displacements  of  938  lbs.  and  1095  lbs.  the  power 
requirements  rise  very  sharply  at  model  speeds  of  3.3  and  3.25  MPH 
respectively,  and  this  sharp  rise  is  correlated  with  the  appearance 
of  a  void  at  the  vehicle  shoulder  and  attendant  poor  flow  to  the 
leading  grousers.  At  a  model  displacement  of  1250  lbs.  (80,000  lbs 
full  scale)  this  phenomenon  occurs  at  a  model  speed  of  3.9  MPH. 

For  a  model  displacement  of  1400  lbs.  the  power  rises  sharply  at 
lower  speeds  for  another  reason,  namely  that  the  bow  of  the  vehicle 
is  completely  covered  with  water  at  3.25  model  MPH.  It  is  curious 
that  for  a  1250  lb.  model  the  onset  of  poor  flow  to  the  grousers 
is  not  reflected  in  curves  of  nominal  propulsive  efficiency  q, 
defined  by 
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Figure  7.1  shows  two  curves  of  q  based  on  HP1  with  tracks 
running  at  zero  slip  and  stationary,  for  1250  lbs.  model  displace¬ 
ment.  Tj  rises  steadily  with  the  tracks  at  zero  slip,  and  rises 
steadily  but  levels  off  at  3  model  MPH  with  tracks  stationary. 

With  938  lb.  model  displacement  we  find  a  sharp  drop  in  nominal 
efficiency  at  a  model  speed  of  3.2  MPH,  which  correlates  well  with 
the  appearance  of  a  shoulder  void  and  a  sharp  rise  in  track  horse¬ 
power  at  3.3  MPH  (see  above). 

We  may  summarize  by  saying  that  on  the  basis  of  data 
used  in  preparing  this  report  neither  of  the  questions  1) ,  2) 

above  can  be  answered  in  a  satisfactory  way. 

The  nominal  efficiency  defined  above  differs  from  the 
efficiency  based  on  the  propulsive  force  exerted  by  the  tracks 
on  the  model.  Reference  [2]  states  that  the  efficiency  (undefined) 
is  of  the  order  of  10-13%  using  the  best  grousers.  This  is  of 
the  same  order  as  what  we  show  in  Figure  7.1  but  the  data  in  [2] 
is  less  variable  with  speed. 

Given  the  present  vehicle  as  a  starting  point  it  is 
clear  from  Figure  7.1,  in  spite  of  the  inaccuracy  of  the  data, 
that  there  is  more  to  be  gained  by  improving  propulsive  efficiency 
than  by  reducing  resistance.  However,  if  propulsive  efficiency 
is  increased  the  relative  importance  of  reducing  resistance  will 
increase.  This  increase  in  the  importance  of  reducing  resistance 
will  be  rapid  since  the  top  speed  of  the  vehicle  will  increase 


IN/  A  SUBSIDIARY  OF  CONTROL  DATA  CORPORATION 


X-100,  TRACKS  STATIONARY, 
LEVEL  TRIM,  938  LBS 


J _ _  1 . 1„  —  ,-.-.1 _ _ I 

.1  .2  .3  .4  .5  1-  V/\/gL 


i 

I 


1— -J _ 1 

2  3  4  MPH  MODEL 


8  MPH  fULL  2CALE 


♦ 


Figure  7.1 

Nominal  Propulsive  Efficiencies  of 
Track- Propelled  X-100  Model,  1250  lbs.  and  938  lbs. 
Model  Displacement,  Level  Trim. 
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for  constant  installed  power,  and  hence  the  drag  will  also 
increase. 

We  believe  that  it  would  be  appropriate  to  continue 
our  development  of  theoretical  models  of  resistance  and  flow. 

At  the  same  time,  the  development  and  test  of  improved 
configurations,  following  the  bulb  whose  calculated  performance 
is  shown  in  Figure  5.1,  should  be  actively  pursued. 
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